Challenges to the DGP Model from Horizon-Scale Growth and Geometry 
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We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self- 
accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, super- 
novae and Hubble constant data by implementing an effective dark energy prescription for modified 
gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between 
distance measures and horizon scale growth in this model. Growth alterations due to perturbations 
propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat 
cosmology, the maximum likelihood DGP model is nominally a 5.3cr poorer fit than ACDM. Cur- 
vature can reduce the tension between distance measures but only at the expense of exacerbating 
the problem with growth leading to a 4.8a result that is dominated by the low multipole CMB 
temperature spectrum. While changing the initial conditions to reduce large scale power can flatten 
the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of 
recent results from WMAP5. The failure of this model highlights the power of combining growth 
and distance measures in cosmology as a test of gravity on the largest scales. 



I. INTRODUCTION 

On the self-accelerating branch of the Dvali- 
Gabadadze-Porrati (DGP) braneworld model [l[ , cosmic 
acceleration arises from a modification to gravity at large 
scales rather than by introducing a form of mysterious 
dark energy with negative pressure Q • In this model our 
universe is a (3 + l)-dimensional brane embedded in an 
infinite Minkowski bulk with differing effective strengths 
of gravity. The relative strengths define a crossover scale 
on the brane beyond which (4 + l)-dimcnsional gravity 
and bulk phenomena become important. 

A number of theoretical and observational problems 
of the self-accelerating branch of DGP have been re- 
cently uncovered. This branch suffers from pathologies 
related to the a ppe arance of ghost degrees of freedom 
d, 0, [H, 0, 0, H, fol, 1 1 0t| . Ghosts can lead to runaway excita- 
tions when coupled to normal modes and their existence 
can invalidate the self-accelerating background solution 
as well as linear perturbations around it. 

Observational problems also arise if one posits that 
ghosts and strong coupling do not invalidate the idea of 
self- acceleration itself, i.e. that our Hubble volume be- 
haves in a manner that is perturbatively close on scales 
near the horizon to the modified Friedmann equation 
specified on this branch. These problems fall into two 
classes, those due to the background expansion history 
and those due to the growth of structure during the ac- 
celeration epoch. 

In a flat spatial geometry the DGP model adds only 
one degree of freedom to fit acceleration, the crossover 
scale. Like ACDM, the extra degree of freedom can be 
phrased as the matter density relative to the critical den- 
sity. Remarkably, with this one parameter, the ACDM 
model can fit three disparate sets of distance measures: 



the local Hubble constant and baryon acoustic oscilla- 
tions, the relative distances to high-rcdshift supernovae 
(SNe), and the acoustic peaks in the cosmic microwave 
background (CMB). The flat DGP model cannot fit these 
observations simultaneously and the significance of the 
discrepancy continues to grow. The addition of spatial 
curvature can help alleviate some, but not all, of this 
tension dUdEl. 

On intermediate scales that encompass measurements 
of the large scale structure of the Universe, the early 
onset of modifications to the expansion also imply a sub- 
stantial reduction in the linear growth of structure which 
is itself slightly reduced by the modification to gravity 
[Til. [To]]. If extended to the mildly nonlinear regime where 
strong coupling effects may alter growth, this reduction 
is in substantial conflict with weak lensing data [l6| . 

On scales approaching the crossover scale, which are 
probed by the CMB, the unique modification to grav- 
ity in the DGP model itself strongly alters the growth 
rate. Here the propagation of perturbations into the bulk 
requires a (4 + l)-dimensional perturbation framework 
[Tn |. The approximate iterative solutions introduced in 
Rfn. [TU have been recently verified to be sufficiently 
accurate by a more direct calculation [l9[ but still are 
computationally too expensive for Monte Carlo explo- 
rations of the DGP parameter space. More recently, a 
(3 + l)-dimensional effective approach dubbed the pa- 
rameterized post-Friedmann (PPF) framework (2^, [2l[ 
has been developed that accurately encapsulates mod- 
ified gravity effects with a closed effective dark energy 
system. The PPF approach enables standard cosmolog- 
ical tools such as an Einstein-Boltzmann linear theory 
solver to be applied to the DGP model. 

In this Paper, we implement the PPF approach to 
DGP and conduct a thorough study of the tension be- 



2 



tween CMB distance, energy density, and growth mea- 
sures, along with the SNe and local distance measures, 
across an extended DGP parameter space. We find that 
even adding epicycles to the DGP model does not sig- 
nificantly improve the agreement with the data. Curva- 
ture, while able to alleviate the problem with distances, 
exacerbates the problem with growth. Changing the ini- 
tial power spectrum to remove excess power in the tem- 
perature spectrum destroys the agreement with recent 
polarization measurements from the five-year Wilkinson 
Microwave Anisotropy Probe (WMAP) [22l |. 

The outline of the paper is as follows. In Sjll] we re- 
view the impact of the DGP modifications on distance 
measures and the growth of structure emphasizing an 
effective dark energy PPF approach that is detailed in 
Appendix A and compared with a growth-geometry split- 
ting approach in Appendix B. We present the results of 
the likelihood analysis in Finl and discuss these results in 

m 



II. DGP DISTANCE AND GROWTH 
PHENOMENOLOGY 

In the DGP model, gravity remains a metric theory 
on the brane in a background that is statistically ho- 
mogeneous and isotropic. Its modifications therefore are 
confined to the field equations that relate the metric to 
the matter. Once the metric is obtained, all the usual 
implications for the propagation of light from distant 
sources and the motion of matter remain unchanged. 
In this section, we review the DGP modifications to 
the background metric, or expansion history, and the 
gravitational potentials, or linear metric perturbations. 
We cast these modifications in the language of an ef- 
fective dark energy contribution under ordinary gravity 



A. Background Evolution 

That the DGP model is a metric theory in a statisti- 
cally homogeneous and isotropic universe imposes a back- 
ground Friedmann-Robertson- Walker (FRW) metric on 
the brane. The FRW metric is specified by two quanti- 
ties, the evolution of the scale factor a and the spatial 
curvature K. The DGP model modifies the field equa- 
tion, i.e. the Fricdmann equation, relating the evolution 
of the scale factor H = a^da/dt to the matter-energy 
content. On the self-accelerating branch it becomes (see, 
e.g., Rfns. UM) 




K 
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where r c is the crossover scale and subscript i labels the 
true matter-energy components of the universe to be dis- 



tinguished below from the effective dark energy contri- 
bution. 

On the other hand, given the metric the matter evolves 
in the same way as in ordinary gravity 
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Overdots represent derivatives with respect to conformal 
time r\ = J dt/a. From these relations, one sees that as 
a — > oo, H — > r" 1 and the Universe enters a de Sitter 
phase of accelerated expansion. 

In the limit that r c — > oo the ordinary Friedmann equa- 
tion is recovered. The effect of a finite r c compared with 
the Hubble scale can be encapsulated in a dimensionless 
parameter f2 rc much like the usual contributions of the 
density fii = SirGpi/ZHQ and curvature Qk = —K/H^ 
to the expansion rate 



(3) 



The modified Friedmann equation {T]) today then be- 
comes the constraint equation 



1 = Q rc +J2^t 
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It is convenient and instructive to recast the impact of 
r c as an effective dark energy component. With the same 
background evolution, the effective dark energy will have 
an energy density of 
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Conservation of its energy-momentum tensor is guaran- 
teed by the Bianchi identities and requires its "equation 
of state" w e = P e / p e to be given by Eq. ^ 



3(H 2 
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If we define f2 e in the same way as fl;, the usual constraint 
condition applies + ^ e + &k = 1- Comparing it to 

Eq. (U)) , we obtain the following relationship between fl e 
and fi r 



fi, 



(7) 



Given w e and f2 e , we can now describe the background 
evolution of the DGP cosmology by using the ordinary 
Friedmann equation for H. Likewise the Hubble param- 
eter specifies the comoving radial distance 



D(z) = I 
Jo 



dz' 
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and the luminosity distance 
(1 + *) 



dhiz) = 
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sin (V^Q^HoD(z) ), (9) 
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FIG. 1: Equation of state of the effective dark energy w e for 
the self-accelerating DGP model with Q m = 0.26 and 51 k = 0. 



as usual. 

The effective dark energy for DGP has quite a dif- 
ferent equation of state from that of a cosmological con- 
stant. For example, when there is no curvature, w e starts 
at — 1/(1 + fi m ) at the present, approaches —1/2 in the 
matter dominated regime and —1/3 in the radiation dom- 
inated regime (see Fig.[TJ. 

With the same values of parameters, DGP will have 
a larger amount of effective dark energy at a given red- 
shift as w e is always larger than —1, and the universe 
will expand at a larger rate than in ACDM. This reduces 
the absolute comoving radial distance to a distant source. 
For relative distance measures like the SNe, the reduction 
in H cIl(z) can be compensated by lowering the frac- 
tional contribution of matter through Q, m . The same is 
not true for absolute distances, such as those measured 
by the CMB, baryon oscillations, and Hubble constant, 
if the overall matter contribution to the expansion rate 
J] m f/nfl -3 remains fixed. 



B. Structure Formation in the Linear Regime 

The same methodology of introducing an effective dark 
energy component for the background expansion history 
applies as well to the linear metric perturbations that 
govern the evolution of large scale structure. 

To define the effective dark energy, one must first pa- 
rameterize the solutions to the (4+ l)-dimcnsional equa- 
tions involving metric perturbations in the brane as well 
as the bulk |17|. Three simplifications aid in this pa- 
rameterization [21| . The first is that at high redshift the 
effect of r c and the extra dimension goes away rapidly [see 
Eq. ([Tjl], The parameterization needs only to be accurate 
between the matter dominated regime and the present. 



Likewise, the predictions of ACDM for the relation- 
ships between the matter and baryon density at recom- 
bination and morphology of the CMB acoustic peaks re- 
main unchanged. As a consequence, the shape of the 
CMB acoustic peaks still constrains the physical cold 
dark matter and baryon density fl c h 2 and fl^h 2 as usual. 
With these as fundamental high redshift parameters, in 
a flat spatial geometry, the DGP degree of freedom r c is 
then specified by either f2 m = 1 — O e = Q c + or Hq. 
We shall see below that in a flat universe the competing 
requirements of CMB and SNe distance measures on r c 
will slightly shift the values of £l c h 2 and tl^h 2 to reach a 
compromise at the expense of the goodness of fit. 

The second simplification is that well below the hori- 
zon, perturbations on the brane are in the quasi-static 
regime where time derivatives can be neglected in com- 
parison with spatial gradients and propagation effects 
into the bulk are negligible. This allows the equations 
to be simply closed on the brane by a modified Pois- 
son equation (|12[) that can be recast as arisin g fr om the 
anisotropic stress of the effective dark energy [l4l [l5j . 

The final simplification is that on scales above the hori- 
zon, the impact of the bulk perturbations on the brane 
becomes scale free and depends only on time fl8j through 
dimcnsionlcss combinations of H , r c and K . Furthermore 
on these scales, generic modifications to gravity are fully 
defined by the Friedmann equation and the anisotropic 
stress of the effective dark energy [29|, H3] ■ 

As shown in Rfn. [2lJ an d detailed in Appendix A, 
interpolation between these two limits leads to a sim- 
ple PPF parameterization of DGP on all linear scales. 
This parameterization has been verified to be accurate 
at a level substantially better than required by cosmic 
variance by a direct computation of bulk perturbations 
[l9l ]. With such a parameterization, efficient Einstein- 
Boltzmann codes such as CAMB [3lJ can be modified to 
calculate the full range of CMB anisotropy, as is done in 
this paper. 

The result of this calculation is that compared with 
ACDM, the growth of structure in the DGP model during 
the acceleration epoch is suppressed. In the quasi-static 
regime, this suppression is mostly due to the higher red- 
shift extent of the acceleration epoch discussed in the 
previous section (see also Fig. [1} with a small compo- 
nent from the effective aniso trop ic stress or modification 
to the Poisson equation [3, [HI- On scales approaching 
r c , the effect of the anisotropic stress becomes much more 
substantial due to the perturbations propagating into the 
bulk [l8| resulting in an even stronger integrated Sachs- 
Wolfe (ISW) effect in the CMB anisotropy power at the 
lowest multipoles [l3j . 



III. CONSTRAINTS FROM CURRENT 
OBSERVATIONS 

In this section, we employ Markov Chain Monte Carlo 
(MCMC) techniques to explore constraints on the DGP 
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parameter space from current observations, and compare 
them with the successful ACDM model. The data we 
use are: the Supernovae Legacy Survey (SNLS) 1331. the 
CMB anisotropy data from the five-year WMAP [33[ for 
both temperature and polarization (TT + EE + TE), 
and the Hubble constant measurement from the Hubble 
Space Telescope (HST) Key Project 0]. 

We use the public MCMC package CosmoMC 
with a modified version of CAMB for DGP described 
in Appendix A, to sample the posterior probability dis- 
tributions of model parameters. The MCMC technique 
employs the Metropolis-Hastings algorithm [36|, [37| for 
the sampling, and the Gelman and Rubin R statistic 
[3^ for the convergence test. We conservatively require 
R — 1 < 0.01 for the eight chains we run for each model, 
and this generally gave us ~ 5000 independent samples. 

The SNe magnitude-redshift relation probes the rela- 
tive luminosity distance between the low and high red- 
shift sample. The luminosity distance itself is completely 
determined by the background expansion history through 
Eq. ©. On the other hand the power spectra of the 
CMB anisotropy probes not only the background expan- 
sion, but also the growth of structure. To better separate 
the two types of information, we also consider a canon- 
ical scalar field ("quintessence" or QCDM) model with 
the same expansion history as DGP. This model is de- 
fined by the equation of state parameter w e and density 
parameter fl e in Eqs. © and (JT]). 

With only scalar perturbations in considera- 
tion, our basic parameter set is chosen to be 
{Clbh 2 , fl c h 2 , 8 S , r, n s , A s }, which in turn stand for 
the density parameters of baryons and cold dark matter, 
angular size of the sound horizon at recombination, 
optical depth from reionization (assumed to be instan- 
taneous), spectra index of the primordial curvature 
fluctuation and its amplitude at fc* = 0.002 Mpc -1 , 
i.e., A 2 C = A^k/k*)"'- 1 . Also, we follow Rfn. and 
include A$z, with flat prior of < A$z < 2, to account 
for the contributions to the CMB power spectra from 
Sunyaev-Zeldovich fluctuations. The lensing effect on 
the CMB is neglected. Note that in all the three models 
we considered, i.e. self-accelerating DGP, ACDM and 
QCDM, we have the same parameter sets and priors 
except that wc restrict the DGP parameter space to 
Hqt c > 1.08 so that metric fluctuations remain well 
behaved [see Eq. I|16p], Wc apply this prior to the 
QCDM model as well for a fair comparison. In practice, 
the excluded models are strongly disfavored by the data 
and the prior is only necessary for numerical reasons. 



A. Flat Models 

We start with the minimal parameterization of a flat 
universe with scale free initial conditions. In this case 
the three model classes ACDM, DGP and QCDM all 
have only one parameter that describes acceleration. In 
the chain parameters this is 9 S but can be equivalently 



TABLE I: Mean and marginalized errors for various parame- 
ters of the self-accelerating DGP, QCDM with the same ex- 
pansion history as the DGP and ACDM models from SNLS 
+ WMAP5 + HST, assuming a flat universe. The first 6 pa- 
rameters are directly varied when running the Markov Chains, 
while the others are derived parameters, as are in the follow- 
ing tables. 



parameters 


DGP 


QCDM 


ACDM 


lOOfith 2 


2.36±0.07 


2.32±0.07 


2.25±0.06 


n c h 2 


0.090±0.005 


0.090±0.005 


0.109±0.005 


1OO0 S 


1.042±0.003 


1.041±0.003 


1.040L0.003 


T 


0.10±0.02 


0.10±0.02 


0.09±0.02 


n s 


1.00±0.02 


0.99±0.02 


0.96±0.01 


In [lO 10 ^] 


3.02±0.05 


3.05±0.04 


3.18±0.04 


Ho 


66 ±2 


65±2 


72±2 




0.26 ±0.02 


0.26±0.02 


0.26±0.02 


n re 


0.136 ±0.009 







TABLE II: Parameters and the likelihood values at the best- 
fit point of the self-accelerating DGP, QCDM with the same 
expansion history as the DGP and ACDM models fitting to 
SNLS + WMAP5 + HST, assuming a flat universe. 



parameters 


DGP 


QCDM 


ACDM 


100n b h 2 


2.37 


2.32 


2.26 


Q.J1 2 


0.0888 


0.0907 


0.110 


1006» s 


1.04 


1.04 


1.04 


T 


0.0954 


0.0998 


0.0825 


n s 


0.998 


0.983 


0.959 


In [10 10 A.] 


3.01 


3.06 


3.18 


Ho 


66.0 


65.1 


71.6 




0.258 


0.269 


0.258 




0.138 






-2hxL 


2805.8 


2797.6 


2777.8 



defined as the derived parameters Hq or f2 m . The con- 
straints on the three model classes are given in Table H] 
for the means and marginalized errors on various param- 
eters. Table |H] shows the best-fit values of the parame- 
ters and the corresponding likelihoods, which serve as a 
"goodness of fit" criterion. 

First, we compare the constraints on the QCDM model 
with those on ACDM. The differences are expected to 
reflect those between the background expansion histories 
of the DGP and ACDM models. In spite of the clustering 
effects on the largest scales of a quintessence dark energy, 
the difference between QCDM and ACDM is completely 
encoded in the equation of state of their dark energy 
components given the fixed sound speed of quintessence. 

Constraints from the SNe magnitudes come from the 
dimcnsionlcss luminosity distance H dz J (z) [see Eq. 
once the unknown absolute magnitude is marginalized. 
In order to match the predictions for HqiIl(z) of a flat 
ACDM model, we would expect that the QCDM model 
has a smaller f2 m to compensate the larger w e its dark 
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TABLE III: Mean and marginalized errors for various param- 
eters of the self-accelerating DGP, QCDM with the same ex- 
pansion history as the DGP and ACDM models from SNLS 
+ WMAP5 + HST, allowing curvature. 



parameters 


DGP 


QCDM 


ACDM 




2.37±0.07 


2.34±0.07 


2.25±0.06 


Q c h 2 


0.096±0.006 


0.098±0.006 


0.108±0.006 


1OO0 S 


1.043±0.003 


1.042±0.003 


1.040±0.003 


T 


0.09±0.02 


0.09±0.02 


0.09±0.02 


«K 


0.019±0.008 


0.027±0.008 


-0.004±0.009 


n s 


1.00±0.02 


0.99±0.02 


0.96±0.01 


In [10 10 A S ] 


3.02±0.05 


3.06±0.04 


3.18±0.04 


Ho 


74±4 


77±5 


70±4 




0.22±0.03 


0.20±0.03 


0.27±0.03 


n rc 


0.15±0.01 







FIG. 2: Predictions for the power spectra of the CMB tem- 
perature anisotropies Cj T of the best-fit DGP (solid), QCDM 
with the same expansion history as DGP (short-dashed), and 
ACDM (dashed, coincident with QCDM at low £) models ob- 
tained by fitting to SNLS + WMAP5 (both temperature and 
polarization) + HST, assuming a flat universe. Bands repre- 
sent the 68% and 95% cosmic variance regions for the DGP 
model. Points represent WMAP5 measurements; note that 
noise dominates over cosmic variance for £ > 500. 



energy has. However, lowering O m also shortens more 
of the distance to the last scattering surface and hence 
increases the angular size of the sound horizon (see mi Ap . 

The physical scale of the sound horizon can partially 
compensate and is controlled by Q, c h? and fl^h 2 but 
these parameters are also well measured by the shape 
of the peaks and can only be slightly adjusted at a cost 
to the goodness of fit. Thus the parameter ranges in 
Table U for fl c h 2 decrease and tti,h 2 increase slightly 
which both have the effect of decreasing the angular 
size of the horizon while £l m remains nearly unchanged. 
This compromise between the energy density and dis- 
tance constraints results in tension between the CMB 
and SNe data. This tension shows up as a differ- 
ence between the — 21ni values of these two models: 
2 In L{ ACDM) - 21nL(QCDM) ~ 20. 

The QCDM model also favors a larger n s and a smaller 
A s . This is a consequence of the larger IS W effect in the 
QCDM model due to a slower growth rate. Tilting the 
spectrum can compensate for the excess power in the 
low-^ modes (as shown by the near-coincidence of the 
short-dashed and dashed curves in Fig. [2]). 

Next we compare the constraints on the DGP model 
with those on QCDM. Since the two models have 
the same expansion history, the differences are entirely 
caused by the differing growth rates. Due to the prop- 
agation of perturbations into the bulk for scales near 
r c and the opposite effect of dark energy clustering in 
QCDM, there is a substantially stronger ISW effect in 



the first few multipoles of the CMB anisotropy power 
[l3j . Since this effect is only important on the largest 
angular modes in the CMB TT power spectra which is 
further limited by the large cosmic variance, the param- 
eter ranges for these two models do not differ signifi- 
cantly. Nonetheless from Fig. [2j it is clear that the best- 
fit DGP model over predicts the low--^ modes anisotropy, 
though as before, n s and A s adjustments try to re- 
duce the primordial perturbations on large scales. This 
leads to DGP being an even worse fit than QCDM with 
21nL(QCDM) - 21nL(DGP) ~ 8. 

When DGP is compared to ACDM, this brings the 
change in — 2 In L for the maximum likelihood parameters 
to ~ 28, where ~ 70% is driven by the background expan- 
sion, while ~ 30% by the dynamical effects on structure 
growth. 



B. Adding in Curvature 

From our analysis in the above section, the flat DGP 
model is a poor fit to the current observations mostly 
because it cannot simultaneously satisfy the geometrical 
requirements of the relative luminosity distances of the 
SNe and the angular size of the sound horizon at recom- 
bination with a single parameter. Since curvature has 
more of an effect on high redshift distance measures, the 
tension in the distance measures can be alleviated by in- 
cluding VLk in the parameter space [pll . [P3T ] . Our results 
arc given in Table [In] and Table [TV] 

With curvature, the — 2h\L of the maximum likeli- 
hood ACDM model almost has no improvement. This 
is consistent with the results of Rfn. [3{|, who found 
strong limits on curvature in ACDM, by fitting WMAP5 
data combined with SNe or HST. As expected the max- 
imum likelihood model in the QCDM space improves in 
— 21nL by ~ 10. The QCDM model needs an open uni- 
verse to increase the distance to last scattering to com- 
pensate the smaller f2 m , consistent with the findings by 
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TABLE IV: Parameters and likelihood values at the best- 
fit point of the self-accelerating DGP, QCDM with the same 
expansion history as the DGP and ACDM models fitting to 
SNLS + WMAP5 + HST, allowing curvature. 



parameters 


DGP 


QCDM 


ACDM 


WOtlth 2 


2.38 


2.36 


2.27 


n c h 2 


0.0937 


0.0960 


0.107 


1OO0 S 


1.04 


1.04 


1.04 


T 


0.0887 


0.0914 


0.0884 


fix 


0.0189 


0.0268 


-0.00553 


n s 


0.996 


0.992 


0.959 


In [W 10 A S ] 


3.02 


3.05 


3.18 


Ho 


73.8 


78.3 


69.8 




0.216 


0.195 


0.266 


fir- a 


0.149 






-2hlL 


2800.8 


2787.2 


2777.5 



[12| . Even with this additional freedom, distance mea- 
sures remain in slight tension due to the Hubble con- 
stant since lowering fl m with fl m h 2 well determined by 
the CMB implies a higher Hubble constant [l3| ■ The al- 
lowed amount of this shift is limited by the HST Key 
Project constraint of Ho = 72 ± 8 km s _1 Mpc -1 . For 
QCDM the baryon acoustic oscillation constraint would 
already disfavor such a shift [4(1 S3, H3, |43[ but its appli- 
cation to DGP requires cosmological simulations of the 
strong coupling regime. Note that with the distance ten- 
sion partially removed, the shifts in Q c h 2 and fl^h 2 are 
reduced. 

The lowering of fi m , in addition to the large equation 
of state parameter of the quintessence, also causes matter 
domination to terminate at an earlier rcdshift, and loads 
to a stronger ISW effect in the QCDM model, as can be 
seen in Fig. [3J again with a partial compensation from 
n s and A s . The net difference of — 21nL values of the 
QCDM model compared to ACDM is ~ 10, 50% smaller 
than before but still a significantly poorer fit. 

The situation is even worse for DGP. Here the enhance- 
ment of the ISW effect at low Vl m is even more substan- 
tial. Thus the mean value of flu is smaller than the 
optimal one for the distance constraints in QCDM. Even 
adjusting the other parameters to give the maximum like- 
lihood model shown in Fig. [3l the poor fit is noticeable 
at the low multipoles. For example the probability of ob- 
taining a quadrupole as extreme as the observations from 
the DGP maximum likelihood model is ~ 1% compared 
with ~ 6% for the ACDM maximum likelihood model. 
The net difference in — 21nL by including curvature as 
a parameter is only ~ 5 in DGP showing only marginal 
evidence for curvature in the model at best. 

Moreover, the difference from ACDM remains substan- 
tial with —2 A In L ~ 23. Since the difference of — 2A In L 
from QCDM can be attributed to the ISW effect, ~ 40% 
of this difference is driven by the background expansion, 
and ~ 60% by the dynamical effects on structure growth. 
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FIG. 3: Predictions for the power spectra of the CMB tem- 
perature anisotropies Cj T of the best-fit DGP (solid), QCDM 
with the same expansion history as DGP (short-dashed), 
and ACDM (dashed) models obtained by fitting to SNLS + 
WMAP5 (both temperature and polarization) + HST, allow- 
ing curvature. 



C. Changing the Initial Power 

The adjustments of A s and n s in the above examples 
suggest that perhaps a more radical change in the initial 
power spectrum can bring DGP back in agreement with 
the data. For example, one can sharply reduce large scale 
power in the temperature spectrum by cutting off the ini- 
tial power spectrum on large scales, below a wavenumber 
&min- While this is a radical modification and has no par- 
ticular physical motivation, it is useful to check whether 
such a loss of power could satisfy the joint temperature 
and polarization constraints. 

CMB polarization at these scales arises from reioniza- 
tion which occurs at a substantially higher rcdshift than 
DGP modifications affect. A finite polarization requires 
not only ionization but also large scale anisotropy at this 
epoch. Eliminating initial power on these scales elimi- 
nates the polarization as well. The EE power in the low 
multipoles has now been measured at 4-5cr level J22| lead- 
ing to a significant discrepancy if the large scale power is 
removed in the model. 

Given that including this parameter does not improve 
the fit, instead of adding it to the MCMC parameter 
space, we illustrate its effects on the maximum likelihood 
DGP model in pIIBl 

By maximizing the likelihood for this model to the 
TT power alone, we find the best agreement is obtained 
at k min = 8 x 10~ 4 Mpc -1 , with -2AlnL TT ~ -12 
compared to fc m ; n = 0. Model predictions with this fc m i n 
are plotted in Fig. [4] together with the WMAP 5 year 
data. 

With this power truncation the over prediction prob- 
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FIG. 4: Predictions for the power spectra of the CMB tem- 
perature anisotropies Cj T of the best-fit DGP model as found 
in Bl without cutting off any large scale primordial per- 
turbations (solid) and with a cut-off scale of fcmin = 8 x 10 -4 
Mpc -1 (dotted) - the best-fit scale obtained when fitting to 
the WMAP 5 year TT data alone, while all other parameters 
are fixed at their best-fit values with fcmin = 0. 



lem for the low-f TT power is alleviated, but the EE po- 
larization power on large angular scales is significantly 
reduced (see Fig. [5]). When the polarization data of 
TE + EE arc included, we find -2AlnL a11 - 6 at 
fcmin = 8 x 10 -4 Mpc -1 , and when fc m i n is varied, the 
combined data does not favor any positive values of fc m ; n 
as shown in Fig. [5] We conclude that though the DGP 
model can be made a better fit to the TT power spectrum 
with a large scale cut-off, polarization measurements are 
now sufficiently strong to rule out this possibility. While 
we have only included instantaneous rcionization mod- 
els, changing the ionization history to have an extended 
high redshift tail can only exacerbate this problem by 
increasing EE power in the i ~ 10 — 30 regime [H|. 



IV. DISCUSSION 

We have conducted a thorough Markov Chain Monte 
Carlo likelihood study of the parameter space available 
to the DGP self-accelerating braneworld scenario given 
CMB, SNe and Hubble constant data. To carry out this 
study, we have introduced techniques for characterizing 
modified gravity and non-canonical dark energy candi- 
dates with the public Einstcin-Boltzmann code CAMB 
that are of interest beyond the DGP calculations them- 
selves. 

We find no way to alleviate substantially the tension 
between distance measures and the growth of horizon 
scale fluctuations that impact the low multipole CMB 
temperature and its relationship to the polarization. In 



FIG. 5: Predictions for the power spectra of the CMB E-mode 
polarization Cf E of the best-fit DGP model as found in mil Bl 
without cutting off any large scale primordial perturbations 
(solid) and with a cut-off scale of fcmin = 8 x 10 -4 Mpc -1 (dot- 
ted) - the best-fit scale obtained when fitting to the WMAP 
5 year TT data alone, while all other parameters are fixed 
their best-fit values. Note here, according to Rfn. [2^], the 
reionization feature at the lowest-^ modes is preferred by the 
data through A\ 2 = 19.6. 




FIG. 6: The total log-likelihood of TT + TE + EE as a func- 
tion of the cut-off scale fc m i n of the primordial perturbations, 
shown as the difference from its value at fc m i n = 0. We find 
that when polarization is included, the combined data does 
not favor any positive values of fcmin- The local minimum 
shows up at the scale of fc m i n = 8 x 10~ 4 Mpc -1 , which is the 
one favored by the TT data alone. Note all other parameters 
are fixed at their best-fit values of the DGP model as found 

in gnrn 



particular we show that the maximum likelihood flat 
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DGP model is a poorer fit than ACDM by 2A InL = 28, 
nominally a (2AlnL) 1 / 2 = 5.3a result. Interestingly, a 
substantial (~ 30%) contribution comes from the change 
in the growth near the crossover scale where perturba- 
tions leak into the bulk. 

Adding in spatial curvature to the model can bring 
the distance measures back in agreement with the data 
but only at the cost of exacerbating the problem with 
growth. The net result is that the maximum likelihood 
only improves by 2A In L = 5 with one extra parameter 
and the difference from ACDM is 2A In L = 23, which is 
still 4.8<7 discrepant with most of the difference arising 
from the changes to growth. 

Furthermore, while the excess power at large angles 
can be reduced by changing the initial power spectrum 
to eliminate large scale power, existing WMAP5 CMB 
polarization measurements already forbid this possibil- 
ity. Specifically, by introducing any finite cut off to the 
initial power spectrum to flatten the temperature power 
spectrum, the global likelihood decreases. 

While it is still possible that the resolution of the ghost 
and strong coupling issues of the theory can alter these 
consequences, it is difficult to see how they can do so 
without altering the very mechanism that makes it a can- 
didate for acceleration without dark energy. The failure 



of this model highlights the power of combining growth 
and distance measures in cosmology as a test of gravity 
on the largest scales. 
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APPENDIX A 

In Appendix A, wc give details of the modification to CAMB employed in the main text to calculate the modified 
growth of perturbations in the DGP model under the PPF prescription. For perturbations of both the metric and 
matter in the various gauges, we use the same notations as in Rfn. [20j. 



1. PPF Description of Modified Gravity in the Linear Regime 



Structure formation in the linear regime for a certain class of modified gravity theories can be equivalently described 
as an effective dark energy component under ordinary gravity. The requirements of this class are that members remain 
metric theories in a statistically homogeneous and isotropic universe where energy-momentum is conserved. This class 
includes the self-accelerating branch of the DGP model. 

Following Rfns. [13, HH, scalar perturbations for the DGP self-accelerating scenario can be parameterized by three 
free functions, g(r],k), f((r]), and /g(^), and one free parameter cr- We shall see in J31 that this is equivalent to 
specifying relationships between the density perturbation, velocity and anisotropic stress of the effective dark energy 
that close the conservation laws of the effective dark energy. 

The metric ratio or anisotropic stress parameter g(rj, k) is defined as 



$+ = g(r), - j^-P T Il T , (10) 

where $ + = ($ + *)/2, $_ = (4> — 1 3/)/2, with $ = 5gij/2gij, ^ = 5 goo / 2#oo as the space-space and time-time 
metric perturbations in the Newtonian gauge, kn = (k/aH), II stands for the anisotropic stress, and the subscript 
"T" denotes sum over all the true components. When anisotropic stresses of the true components are negligible, g 
just parameterizes deviations from GR in the metric ratio of $ + to 

There are two key features of the PPF parameterization and these define the two additional functions f((j]) and 
fciv)- The first is that the curvature perturbation in the total matter comoving gauge £ is conserved up to order k\ 
in the super-horizon (SH) regime in the absence of non-adiabatic fluctuations and background curvature 

t AP T - \c k Pt^t K 1 
lim ^- = 3 K - T^k H V T + -c K fdv)k H V T , 11 

where (77) parameterizes the relationship between the metric and the matter. The second is that in the quasi-static 
(QS) regime, one recovers a modified Poisson equation with a potentially time dependent effective Newton constant 



AttG AtPt + ckPt^It 
fcff>i 1 ~ " c K k 2 H H 2 1 + f G {ri) 



= — una ■ (12) 



Note that even if fa = 0, the Poisson equation for ^> may also be modified by a non-zero g(rj, k). In the above two 
equations, A, V and AP (7^ PA) are density, velocity, pressure perturbations in the total matter comoving gauge, 
and Cjc = 1 — 'SK/k 2 . To bridge the two regimes, an intermediate quantity T is introduced, 

+ r = 4?F 2 G (A tP t + c K P T U T ) , (13) 
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and an interpolating equation is adopted to make sure it dynamically recovers the behavior specified by Eqs. (fTTj) and 



(1 



c T k H ) 



r 



s. 



(14) 



Here the free parameter cr gives the transition scale between the two regimes in terms of the Hubble scale, and the 
source term S is given by, 



s= 9/(aH)-2 9 q> 
.9 + 1 



(g + l)k%H* 



9 



(PtTIt) 
aH 



P T U T 



[(9 + k + 9k)(PT + Pt) - ( P e + Pe)] k H V T . (15) 



By interpolating between two exact behaviors specified by functions of time alone f({rj) and Jg(v)i the PPF param- 
eterization is both simple and general. Moreover, the one function of time and scale 5(77, k) also interpolates between 
two well defined functions of time alone for many models, including DGP. The same is not true of parameterization^ 
that involve the effective anisotropic stress alone or the metric functions directly. 



PPF Parameterization for DGP 



The PPF parameterization for the self-accelerating DGP is given in Rfn. [2l|, which we summarize as the follows. 
On super-horizon scales, the iterative scaling solution developed in Rfn. (l8| is well described by 



9sh(v) = 



9 



1 



0.51 



8Hr c - 1 V Hr c - 1.08, 

and f((ri) = OAgsniy), while in the QS regime, the solution is parameterized by [if 

-1 

H \ 

1 - 2Hr r 1 



gQsiv) = - 3 



3aH 2 



and fciv) ~ 0- O n an arbitrary scale in the linear regime, g is then interpolated by 



9(il, k) 



ffSH + 9Qs(c g k H ) n " 
1 + (cgkn)^ 



(16) 



(17) 



(18) 



with c g = 0.14 and n g 



3. This fitting formula has been shown to give an accurate prediction for the evolution of <f>_ 
according to the dynamical scaling solution [21j . The transition scale for T is set to be cr = 1, i.e. at the horizon scale. 
We note here that the solutions developed in Rfns. [TU, [l8| are for flat universes, so strictly speaking, the above PPF 
parameterization only works for DGP with Q k 0. However for the small curvatures that are allowed by the data, 
we would expect its effect on structure formation to be small and arise from terms such as H 2 — > H 2 + K/a 2 (see, 
e.g., Rfn. [45j). Given the cosmic variance of the low-£ multipoles, these corrections should have negligible impact on 
the results. 



3. "Dark Energy" Representation of PPF 



By comparing the equations that the PPF quantities satisfy with their counterparts in a dark energy system under 
general relativity, we obtain the following relations for the perturbations of PPF's corresponding effective dark energy. 
These relations act as the closure conditions for the stress energy conservation equations of the effective dark energy. 
The first closure condition is a relationship between the PPF T variable and the components of the stress energy 
tensor of the effective dark energy 

Pe A e + 3(p e + p e )^^ + c K p e n e = -^^r. (19) 

ku AirGa^ 
The second closure condition is a relationship for the anisotropic stress 

^n e = -%^$-. (20) 
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Stress energy conservation then defines the velocity perturbation 



V e -Vr 
k H 



with 



H 2 



9 



^G{p e + P e ) F 



S-T 



r 



k 



4irG(p 7 



H 2 



1 



UirGa 2 
k 2 c K 



{g + l)(p T + Pt) ■ 



P T ) V T 



(21) 



(22) 



For details of all the derivations, sec Rfn. [20(. The PPF equation for V then replaces the continuity and Navier-Stokes 
equations for the dark energy. Note that the effective dark energy pressure perturbation, which obeys complicated 
dynamics to enforce the large and small scale behavior, is not used. 



4. Modifying CAMB to Include PPF 

Given the dark energy representation of PPF, we only need to modify the parts in CAMB where dark energy 
perturbations appear explicitly or are needed, in addition to its equation of state which will be specified by the 
desired background expansion. These include the Einstein equations and the source term for the CMB temperature 
anisotropy. Since CAMB adopts the synchronous gauge, in this section, we will express everything in the same gauge. 

The two Einstein equations used in CAMB are, 



^ " c K k„7 lT 



4ttG 
k H H 2 



{ PT 5 S T + Pe S s e ) , 



krjT - K [kL+ 2 ^ T) = AirGa 2 [(p T + P T )v s T + (p e + P e )v s e 



(23) 



(24) 



where h^, rjr are metric perturbations, 5, v are density and velocity perturbations, and superscript "s" labels the 
synchronous gauge. Here, we only need to provide 8 s e and v s e . Given the gauge transformation relation for velocity 



V = v s 



with a = (fiL + 67/r)/2fc 2 , the following expression for v s e is easily obtained from Eq. (j2"T|) 



(p e +P e )v s e = ( Pe + P e )v s T - 



k H H 2 (g + 1) 
4ttG F 



k 



AttG(pt + Pt){v s t + ka) 
kUT 2 



(25) 



(26) 



In order to calculate u|, we need to evaluate a, which we find, with the help of the first closure condition and the two 
Einstein equations, to be given by 



ka = ktfTjT 



AttG 
CKknH 2 



p T 8j 



3(pt + Pt)v s t 



ku 



AttG 



2^n e - k H T 



(27) 



Here, the anisotropic stress is gauge- independent, and is given by Eq. (|20[) . where <&_ is given by gauge-transforming 
the density perturbation in Eq. (|13p according to 



pA = pS s - p- 



(28) 



In addition to a, we need also to specify S and f in order to get v*. Given 5>_, S can be calculated by gauge- 
transforming Vt in Eq. (TTS")) . and T then follows from Eq. P^|) . Provided v s e and P e n e , 8 s e can be obtained by 
gauge-transforming the first closure condition 



p e 5 s e = -c K P e IL e - 3(p e + P e 



c K k 2 H H 2 



kn inG 

The source term for the CMB temperature anisotropy is given by [46| 



(29) 



s T (v,k) = g[ A 



TO 



2a 



E 3E 
46 + 4fc2& 



+ e _K (ryT + a) + Q\ -± + a + 



3S 
2k 2 ! 



3gS 
Ak 2 b ' 



(30) 
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FIG. 7: A comparison of the PPF prediction (upper solid line) with the prediction by parameter-splitting (lower dashed line) 
for the DGP model. 



where the visibility function Q = — /icxp(— k), with k(t]) the optical depth from today to 77 for Thomson scattering, 
At( is the multipole of the temperature anisotropy, is the baryon velocity, £ = At2 — I22AP2, with 2AP2 the 



quadrupole of the polarization anisotropy (for more explicit definitions for Ati and 2Apg, see Rfn. |46|), and b 2 = ck- 
Here we need to specify a and a. a is given by the Einstein equation 



a - 



2— a — 7]t = 



8nG 



(P T n T + P e Tl e ) 



(31) 



Differentiating this equation with respect to 77 also gives us a, which will need the derivative of the effective dark 
energy's anisotropic stress. From energy-momentum conservation and the first closure condition, we obtain 



(P e n e 



h L (Pe+Pe) 




6 +k 3 f h + H 

kn k \ a H 



(32) 



which completes our modification of CAMB. 



APPENDIX B 



In Appendix B, we briefly contrast the PPF prediction for the CMB temperature power spectrum with an attempt to 
approximate the DGP prediction via the parameter-splitting technique (la. l47j. The technique works by splitting the 
dark energy equation of state w into two separate parameters, with one, u> g0O mctry, determining geometric distances 
and the other, ui gr0 wth, determining the growth of structure. We choose u> ge0 mctry = —0.7 and tfgrowth = —0.57, 
which are obtained respectively by fitting to H{z) (for z = to 10) and fitting to the quasi-static growth factor 
(for z = to 2) according to the best-fit flat DGP model (Table HI)) . The result is shown in Fig. [71 One can see 
that parameter-splitting falls short of the PPF prediction on large angular scales. This illustrates the importance 
of correctly modeling the perturbation growth on horizon scales. The parameter-split of w does not contain enough 
freedom to describe perturbation growth in both the sub-horizon (quasi-static) and super-horizon regimes. 



